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ABSTRACT 

V 393 Scorpii is a Double Periodic Variable characterized by a relatively stable non- 
orbital photometric cycle of 253 days. Mennickent et al. argue for the presence of a massive 
optically thick disc around the more massive B-type component and describe the evolutionary 
stage of the system. In this paper we analyze the behavior of the main spectroscopic opti- 
cal lines during the long non-orbital photometric cycle. We study the radial velocity of the 
donor determining their orbital elements and find a small but significant orbital eccentricity 
(e = 0.04). The donor spectral features are modeled and removed from the spectrum at ev- 
ery observing epoch using the light-curve model given by Mennickent et al. We find that the 
line emission is larger during eclipses and mostly comes from a bipolar wind. We find that 
the long cycle is explained in terms of a modulation of the wind strength; the wind has a 
larger line and continuum emissivity on the high state. We report the discovery of highly vari- 
able chromospheric emission in the donor, as revealed by Doppler maps of the emission lines 
Mg II 4481 and CI6588. We discuss notable and some novel spectroscopic features like dis- 
crete absorption components, especially visible at blue-depressed 1 7773 absorption wings 
during the second half-cycle, Balmer double emission with V^/i?-curves showing "Z-type" 
and "S-type" excursions around secondary and main eclipse, respectively, and H/3 emission 
wings extending up to ± 2000 km s _1 . We discuss possible causes for these phenomena and 
for their modulations with the long cycle. 

Key words: stars: early-type, stars: evolution, stars: mass-loss, stars: emission-line, stars: 
winds-outflows, stars: binaries-eclipsing, stars: variables-others 



1 INTRODUCTION 

V393 Scorpii is one of the Galactic Double Periodic Variables 
(DPVs), a group of interacting binaries showing a long photomet- 
ric cycle lasting roughly 33 times the orbital period (Mennickent et 
al. 2003, 2005, 2008, Mennickent & Kolaczkowski 2010, Radek et 
al. 2010). DPVs have been interpreted as semi-detached interact- 
ing binaries ongoing cyclic episodes of mass loss into the interstel- 
lar medium (Mennickent et al. 2008, Mennickent & Kolaczkowski 
2010). The long photometric cycle of V393Sco (253 days) was 
discovered by Pilecki & Szczygiel (2007) by inspecting the ASAS 
catalogue of eclipsing binaries with additional variability. Ultra- 
violet spectra of V393 Sco were studied by Peters (2001) who 
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found evidence for a hot temperature region that should be the ori- 
gin of superionized lines like N V, C IV and Si IV. These lines are 
likely produced by resonance scattering in a plasma of tempera- 
ture T ~ 10 5 K and electron density N e ~ 10 9 cm -3 (Peters & 
Polidan 1984). Broad band photometry was modeled by Mennick- 
ent et al. (2010a, hereafter M10), who found a distance d = 523 
pc and reddening E(B — V) = 0.13. They interpreted observed 
asymmetries in UV lines of highly ionized metals as evidence for a 
high-latitude wind. They also studied orbitally-resolved high reso- 
lution infrared spectra and found a remarkable depressed blue wing 
in the He 1 1083 nm line near secondary eclipse, that interpreted as 
evidence of mass loss through the L3 point. The relatively stable 
ultraviolet spectra during the long cycle yielded these authors to 
suggest that the mass outflow related to the long cycle should be 
through the equatorial plane. 

Posteriorly, Mennickent et al. (2012, hereafter Ml 2) disentan- 
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gled the orbital and long-cycle light curves and performed a multi 
parametric fit to the orbital light curve considering two stellar com- 
ponents plus a circumprimary disc. The best fit is obtained with a 
Roche-lobe filling A7 giant (hereafter the donor or secondary star) 
and a B4 dwarf (the gainer or primary) that is surrounded by an op- 
tically thick disc with a radius almost twice the radius of its central 
star. The disc has 16.600 K in the inner parts and 8.600 K in the 
outer regions. The orbit is circular with separation 35 Rq and the 
inclination is 80 degree. M12 found two bright spots in the disc, that 
are needed to fully reproduce the shape of the orbital light curve. 
These spots are 20% hotter that the disc and face the observer at 
orbital phases 0.90 and 0.45. 

Ml 2 performed a \ 2 minimization procedure to find the sys- 
tem and orbital parameters that best match a grid of theoretical 
evolutionary tracks for binary stars including non-conservative pro- 
cesses. After performing this procedure, Ml 2 find an age of 70 
Myr, but more importantly, they found that the system is just af- 
ter a mass transfer burst that transferred 4 Mq from the donor to 
the gainer in just 400.000 years. The system and orbital parameters 
are reproduced relatively well, but a mismatch between the gainer 
dynamical mass and its luminosity is found. In order to explain this 
discrepancy, M12 argue that not all transferred mass has been accu- 
mulated by the (possibly critically rotating) gainer, but a significant 
fraction (~ 2 Mq) still remains in the disc. Regarding the long cy- 
cle, M 1 2 find that is not produced by changes in the relatively stable 
disc, and invoke a variable high latitude wind as probable cause for 
the phenomenon. Emission in this wind should explain, according 
to them, the non-eclipsing nature of the long variability and their 
red colors at maximum. 

With the aim of studying the nature and origin of the long pho- 
tometric cycle, we carried out a spectroscopic campaign of V 393 
Sco using several high resolution spectrographs at different obser- 
vatories in Chile. Additional spectra taken 20 years ago at Kitt Peak 
National Observatory were also considered. The analysis of the 
main spectral features of V 393 Sco is presented in this paper. We 
postpose for a future paper the analysis of the weaker spectroscopic 
features. 

In this paper long-cycle maximum and minimum are in- 
distinctly named high and low state, respectively. We use the 
ephemerides for the orbital and long cycles given by Ml 2, viz., 
the main orbital photometric minimum happens at HJD — 
2 452 507.7800 + 7.7125772B (Kreiner 2004) and the maximum 
of the long cycle occurs at HJDi = 2452522 + 253E. The ob- 
servations, instruments and reduction techniques are presented in 
Section 2, the results are given in Section 3, a general discussion 
is presented in Section 4 and the conclusions are given in Section 
5. A brief preliminary report of this research was published in a 
conference proceeding by Mennickent et al. (2010b). 



2 SPECTROSCOPIC OBSERVATIONS AND DATA 
REDUCTION 

We conducted spectroscopic observations of V 393 Sco during 2 
years obtaining optical spectra with resolution R ~ 40.000 with 
several echelle spectrographs in Chile. The instruments used were 
CORALIE, HARPS and FEROS (La Silla ESO Observatory, 106, 
6 and 14 spectra, respectively), UVES (Paranal ESO Observatory, 
359 spectra) and DuPont-echelle and MIKE (Las Campanas Ob- 



servatory, 42 and 17 spectra, respectively J] The studied spectral 
regions were 3865-6900 A (CORALIE), 3780-6910 A (HARPS), 
3560-9215 A (FEROS), 3600-9850 A (DuPont-echelle) 3350- 
5000 & 4900-9500 A (MIKE), 3760-4985, 5682-7520 & 7660- 
9464 A (UVES). We also included 9 medium resolution (R « 
15.000) spectra in the Ha region and 2 in the He 4471 region (R « 
22.000) obtained between years 1987-1990 with the COUDE spec- 
trograph in the Kitt Peak National Observatory. Since the long- 
cycle phases are quite uncertain these spectra are only phased with 
the orbital cycle. In general, their variability is not different from 
spectra acquired 20 years later. A total of 555 spectra were acquired 
during our campaign, most of them with signal to noise ratio larger 
than 100. The histograms for the orbital and long-cycle phases of 
the high resolution spectra obtained in the Ha region show that 
our observations cover relatively well the short and long periodic- 
ity (Fig. 1). 

All spectra discussed in this paper are corrected by earth trans- 
lational motion and normalized to the continuum. The RVs are 
heliocentric ones. Reductions were done with IRAlj^] following 
usual procedures for echelle spectrography, including flat and bias 
correction, wavelength calibration and order merging. Occasional 
cleaning of spectra was needed due to the presence of bright fea- 
tures produced by incident cosmic rays or the presence of sharp 
telluric lines. In order to remove these features, we replaced de- 
viant points from a high order fit to the flux distribution in a short 
spectral segment with the flux corresponding to the fitting function. 
Most cosmic ray features and telluric lines were efficiently removed 
with this technique. 

Some of the DuPont spectra showed artifacts around short 
spectral regions produced by interference of light reflected in the 
detector inner layers. These "fringes" were identified and the fluxes 
in the affected regions were interpolated and not considered for the 
analysis. 

The spectra obtained with the CORALIE spectrograph are not 
sky-subtracted. This condition imposed by the instrumental setup 
has not effect for radial velocity and line strength measurements, 
since V 393 Sco is very bright even at full moon and we do not flux- 
calibrate our spectra. Details for our observational runs are given in 
Tables 1, 2 and 3. 



3 RESULTS 

3.1 General spectrum appearance and interstellar reddening 

The optical spectral region of V 393 Sco shows Balmer, He I, Si II 
lines in absorption and sometimes strong emission in Ha, suggest- 
ing a B-type primary surrounded by a variable gaseous envelope, 
along with metallic absorption lines, interpreted as signatures of the 
A-type secondary star. A rapid inspection showed several absorp- 
tion helium lines varying in width and depth, a fact that is consistent 
with a non-photospheric origin. The above description suggests the 
presence of circumstellar matter and it is consistent with the ear- 
lier classification of intermediate-mass interacting binary for this 
system. 

1 Technical descriptions for these spectrographs and their cameras can be 
found in www.lco.cl/ , www.eso.org/sci/facilities/paranal/instruments/ and 
www.eso.org/sci/facilities/lasilla/instruments/. 

2 IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 



A cyclic bipolar wind in the interacting binary V393 Scorpii 3 



Table 1. Summary of LCO and ESO non-UVES spectroscopic observations. N is number of spectra. The HJD at mid-exposure for the first spectrum of the 
series is given. <I> and <!>; refer to the orbital and long-cycle phase, respectively. 



UT-date 


Observatory /Tele scope 


Instrument 


N 


exptime (s) 


HJD-MID (start) 






9007-07-1 Q 

ZUU 1 U 1 - 17 




MIKE 




20 


?4S4^ni 47395 


0.5674 


fl 0335 






MIKE 


5 


230 


7454^01 50990 


5720 


fl 0336 

V..UJJU 


9007 07 9 1 

ZUU 1 <J 1 - Z 1 


T fO/M lapllnn Plnv 


MIKE 


4 


40 


9zlSzl3n9 S3377 
ZH-J'tjUZ.jjj / / 


7048 


0377 


9007-08-0 1 
ZUU 1 -uo-u 1 




MIKE 




920 


94S4313 ft9004 


1422 


0816 


9007 08 09 
ZUU 1 -Uo - UZ 




MIKE 


j 


390 


9zlSzl31zl SS7ft1 


2638 


08 S 9 
U.UoJZ 


9007 08 OS 
ZUU 1 - Uo-UJ 


ij^w/iviagciittii y^itxy 


MIKE 


j 


250 


9zlSzl3 1 7 70700 


6721 


0Q77 
U.UV / / 


9007 08 1 ft 
ZUU / -Uo- 1 


T fO/M iCTpllrin Plnv 


MIKE 


1 


270 


?4Szl398 S8438 

ZtJ^JZoJOtJO 


089 s 

U.UoZJ 


140ft 

U. 1M-UD 


9007 1 1 08 
ZUU i-l i-uo 


T fO/M lapllnn Plnv 
lji_w/lvldgclldll i_ldy 


MIKE 


j 


120 


?4SzL11 9 S07S0 


Qft38 


4724 


9007 1 1 OQ 
ZUU I-l i-uv 


T fO/M iCTpllrin Plnv 
ljL-.W/lvldgclldIl i_ldy 


MIKE 


j 


160 


?4SzL11 3 4Q779 


0Q99 
u.uyzz 


47ft3 


9008 09 94 


T PO/Hn Pnnt 


Fr* n*=M 1 f» 


o 


500 


94S4S90 831 87 


00Q0 


QOOS 
U.VUUJ 


9008 09 9 S 


Ij^W/ULI Jr^UIlL 


Echelle 




500 


?4SzlS91 83S14 

ZtJtJZ l.ojJlt 


1 3Q 1 


Q04S 
U.7U4J 


9008 Oft 1 S 




Echelle 


1 1 


500 


9zlSzlft39 S 1 3ftft 


4895 


3490 
U.J)f zu 


900Q OS 1 4 


T PO/Pln Pnnt 
1j*_W/ LJU Jr^UIlL 


Echelle 




420 


9zlSzlQftS ftSQ37 

ZtJt7UJ.UJ7J / 


ft84ft 


6587 


900Q Oft 1 ft 


T PO/Dn Pnnt 


Echelle 




540-600 


9zlSzlQQ8 ft30 1 Q 
ZH , jH , 77o.OjU 1 y 


9595 


7891 


900Q 08 9S 
ZUU7 uo-ZJ 


T PO/Dn Pnnt 
1j\_W/ U LI r/UIlL 


Echelle 


4 


540 


9zlSS0ft8 S788Q 


09 8Q 

U.UZ07 


OftSS 
U.UUJ J 


900Q 09 9ft 


F^O/ 3 ftm 


HARPS 


j 


300 


94S4888 8771 Q 


79Q1 
U. / -Cy 1 


3SS9 

U.JJJZ 


900Q 09 98 


pen/ 3 ftm 


HARPS 


j 


180 


?4Szl8Q0 89309 
ztJtoyv-.ozjuz 


Q8 1 4 
U.yo It 


3ft9Q 


901 OQ OQ 


F^fV 3 ftm 
HjKJI J.UII1 


HARPS 


j 


180 


94SS-14.8 S1 193 

ZH-JJHtO.J 1 IZj 


9Q00 


Sft79 
U.JO / z 


901 OQ 97 


F^O/ 3 ftm 


HARPS 


j 


o 


?4SS-lft7 473S8 


7480 


6422 


901 OQ 98 
ZU 1 U-UV Zo 


pen/ 3 ftm 


HARPS 


j 
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9zlSS-lft8 48QQ--1 


8 8 Of) 


6462 


901 OQ 30 


pen/ 3 ftm 


HARPS 


j 


250 


94.SSzl.70 47S88 
ZH-J Jt /U.t / JOO 


1380 


6541 


9008 03 1 ft 


pcn/9 9m 


Feros 


'j 
j 


250-300 


9zlSzlS--l1 8Qft07 


7401 


9838 


9008 03 1 7 

ZUUo - UJ- 1 / 


pcn/9 ?m 
CJ\Jl z.zm 


Feros 




300^100 


94S4S--19 Q0141 


8705 


9878 


9008 03 1 8 


pcn/9 ?m 


Feros 




300-400 


?4SzlS--l3 8Q991 


QQ8Q 

U.7707 


9917 


9008 03 1 Q 


F<nO/9 9m 


Feros 


2 


500 


?zlSzlS4zl 8Q83Q 
Z4-JtJt4.070J7 


1 9Q4 

U.1Z74 


QQSft 


9008 03 90 
ZUUO-uj-ZU 


F^O/9 9m 
JZ,iJ vJ/ Z.ZII1 


Feros 




400-500 


?4SzlSzlS 8Q87Q 


9SQ1 

U.ZJ7 1 


QQQft 

U.7770 


9008 Ozl Oft 
zuu o - ut-uu 


F^O/FTTT FR 


rno AT TF 


O 


900 


?4SzlSft9 89SS1 
Z'tJ'I-jOZ.oZjj 1 


4S38 


OftftS 

U.UUOJ 


9008 Ozl 07 


F^O/FTTT PR 

1 L. l W/ 1 L U IjIjK. 


rno AT TF 


o 
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?4SzlSft3 771 OS 
Z4J'tJUJ. / / lOJ 


S7ft4 
U.J 


0709 
u.u /uz 


9008 Ozl 08 


F^O/FTTT FP 


rno AT TF 


4 


900 


?4SzlSft3 771 OS 


S7ft4 

U.J / OH- 


0709 
u.u / uz 


9008 Ozl OQ 


F^O/FTTT FP 
1 lo W/ 1 L U l_,dK 


rnc AT TF 


4 


1000 


?4SzlSftS 831Q0 


843ft 


0784 
u.u / 


9008 Ozl 1 n 


pcn/pTTT pp 


rno AT TF 


j 


900 


94S4.Sftft Q93QS 

Z4-J'1 , J00.7ZJ7J 


9852 


0827 


9008 OS 93 
ZUUo UJ Zj 


F^O/FTTT FP 


rno AT TF 


2 


1000 


94SzlftOQ RAA1A 


5502 


2524 


9008 OS 94 
ZUUo-UJ-Zt 


F^O/FTTT FP 


rno AT TF 

W Ixrt 1j1 C 


4 


1000 


9zlSzlft 1 78498 


6721 


2561 


9008 OS 9S 

ZUUo-UJ-ZJ 


F^O/FTTT FP 

CJU/ Cj U IjIjK 


rnc AT TF 

L V / l\ rt 1 v 1 1 L 


j 


1000 


?4S4ft 1 1 808Q7 


804Q 


9ft01 

u.zoui 


9008 OS 9ft 


F^O/FTTT FP 

1 L. i W/ 1 L U IjEJx 


rno AT TF 


7 


1000 


?4S4ft 19 77847 


03 Oft 

U.7JUU 


9ft3Q 
u.zo jy 


9008 08 1 Q 


F^O/FTTT FP 


rno AT TF 

l_.WK.rt.ljlC 


J 


600 


?4S4ftQ8 S3Sft8 


04Q8 

U.Uf 70 


ftf>9Q 
U.DUZ7 


9008-08-90 
ZUU o uo zu 


F SO/FT TT FR 

LiJL// 1 ,1)1 ,| .IX 


TOR AT TF 


J 


600 


94S4ftQQ 48Qft4 


0.1734 


0.6067 


9008 08 91 


F^O/FTTT FP 
1 L. i W/ 1 1 U 1 j 1 1 l^v 


rno AT TF 

l_.WKrtl_jl.Cj 


J 


600 


94S4700 478Q7 


301 7 
U. JU 1 / 


ft 10ft 
U.DIUD 


9008-08-99 

ZWO UO zz 


F SO/FT TT FR 


TOR AT TF 


2 


600 


94S470I ft3949 

_.7J7 / 1/ 1 .UJZtZ 


0.4513 


0.6151 


2008-08-23 


ESO/EULER 


CORALIE 


2 


600 


2454702.64402 


0.5824 


0.6191 


2008-08-24 


ESO/EULER 


CORALIE 


2 


600 


2454703.59256 


0.7054 


0.6229 


2008-08-25 


ESO/EULER 


CORALIE 


1 


600 


2454704.48019 


0.8205 


0.6264 


2008-08-26 


ESO/EULER 


CORALIE 


1 


600 


2454705.50720 


0.9537 


0.6305 


2008-08-28 


ESO/EULER 


CORALIE 


1 


600 


2454707.51773 


0.2144 


0.6384 


2008-10-01 


ESO/EULER 


CORALIE 


6 


600-800 


2454741.53378 


0.6248 


0.7729 


2008-10-02 


ESO/EULER 


CORALIE 


3 


900-1000 


2454742.52995 


0.7540 


0.7768 


2008-10-03 


ESO/EULER 


CORALIE 


4 


900 


2454743.52755 


0.8833 


0.7807 


2009-04-15 


ESO/EULER 


CORALIE 


1 


1800 


2454936.90276 


0.9560 


0.5451 


2009-04-16 


ESO/EULER 


CORALIE 


3 


1300 


2454937.87533 


0.0821 


0.5489 


2009-04-17 


ESO/EULER 


CORALIE 


4 


1000 


2454938.87689 


0.2120 


0.5529 


2009-05-17 


ESO/EULER 


CORALIE 


11 


900 


2454968.61926 


0.0683 


0.6704 


2009-05-18 


ESO/EULER 


CORALIE 


9 


900 


2454969.72474 


0.2117 


0.6748 


2009-05-19 


ESO/EULER 


CORALIE 


12 


800 


2454970.59597 


0.3246 


0.6782 



The interstellar contribution to the color excess can be deter- 
mined through the analysis of diffuse interstellar bands (DIBs; Mu- 
nari 2000, Weselak et al. 2008). We measured equivalent widths of 
DIBs located at 5780 A, 5797 A and 8620 A, estimating E(B-V) 
= 0.15 ± 0.05, in agreement with the value E(B-V)= 0.13 ± 0.02 
obtained by M10 from the fit to the spectral energy distribution. The 
above probably excludes significant amounts of dust in the line of 



sight to V 393 Scorpii and suggests that the color excess is mainly 
due to interstellar reddening rather than circumstellar matter. 

3.2 The donor's spectrum and its stellar parameters 

In order to find the donor spectral type we compared the spectrum 
observed at minimum light with a grid of synthetic spectra and 
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Table 2. Summary of KPNO spectroscopic observations obtained with the 
Coude spectrograph. One spectrum was obtained per night. <E> is the orbital 
phase. 



Table 3. Summary of observations obtained with the ESO UVES between 
April 02 and October 01, 2008. N is the number of spectra. The exposure 
time was 35 s. The complete table is given in electronic form. 



UT-date 


A A (A) 


exptime (s) 


HJD 


3>o 


1987-05-03 


6520-6695 


1800 


2446918.93455 


0.3595 


1987-08-24 


6520-6695 


2700 


2447031.69913 


0.9804 


1989-04-20 


6520-6695 


2400 


2447636.97037 


0.4589 


1989-04-21 


6520-6695 


2700 


2447637.92552 


0.5827 


1989-04-23 


6520-6695 


2700 


2447639.97292 


0.8482 


1989-04-24 


6520-6695 


3300 


2447640.94252 


0.9739 


1990-08-19 


6520-6695 


2000 


2448122.68538 


0.4359 


1990-08-20 


4420-4535 


2700 


2448123.65323 


0.5614 


1990-08-21 


4420-4535 


3600 


2448124.70935 


0.6983 


1990-08-22 


6520-6695 


1800 


2448125.65894 


0.8214 


1989-94-21 


6520-6695 


2700 


2447637.92552 


0.5827 



looked for the best match in a region dominated by donor spec- 
tral features. For the reference spectrum we used the FEROS spec- 
trum at <fr = 0.000 obtained on March 18, 2008 (hereafter $ D and 
<bi stand for orbital and long cycle phase respectively). The region 
used for the fit, deployed of hydrogen and helium lines but with 
several metallic lines, was the region between 4150 and 4290 A. 

To determine the grid of synthetic fluxes we used atmospheric 
models computed with the line-blanketed LTE ATLAS9 code (Ku- 
rucz 1993), which treats line opacity with opacity distribution func- 
tions (ODFs). The Kurucz's models are constructed with the as- 
sumptions of plane-parallel geometry and hydrostatic and radiative 
equilibrium of the gas. The synthetic spectra were computed with 
the SYNTHE code (Kurucz 1993). Both codes, ATLAS9 and SYN- 
THE were ported under GNU Linux by Sbordone et al. ( 2005} and 
are avai lable onlin^B The atomic data were taken from Castelli & 
Hubrig (20040 The theoretical models were obtained for effective 
temperatures from 6000 to 10000 K with steps of 100 K and for 
surface gravities from 2.0 to 4.5 dex with the step of 0.1 dex. So- 
lar and 0.5 dex higher metallicities were taken into account. The 
grid of synthetic spectra was calculated for five different rotation 
velocities, virSini = 0, 25, 50, 75 and 100 km s _1 . 

We subtracted the template from every grid spectrum, allow- 
ing small wavelength adjustments and analyzed the residuals of the 
resulting difference spectra. The model with the lower root mean 
square had T2 = 7900 K, \0gg2 = 3.0, virsini = 75 km s _1 
and solar metallicity. Several solar metallicity models with temper- 
atures in the range 7500-8100 K showed similarly small residu- 
als. Some higher metallicity models also fit the spectrum but with 
higher temperatures, the best one with 8400 K, logg = 3.2, and 
virSini = 75 km s _1 . We preferred the lower temperature model 
for consistency with our earlier calculation and previously pub- 
lished results (M10 and references therein). 

From the above we estimate T2 = 7900 t.400 K. Our low tem- 
perature limit is consistent with the absence of the G-band in our 
spectra. The low sensitivity of spectral features to the donor temper- 
ature in the chosen spectral region precludes a better determination 
of T2. Other potential temperature spectral diagnostics like calcium 
H & K lines could not be used due to contamination with features 
of other system components. 

A refinement of v-irsini was done by interpolating the grid of 
models to the FWHM of lines observed on main eclipse, obtain- 



N 


AA 
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(A) 


(2450000 +) 


119 


3760-4985 


4558.90352-4740.54141 


120 


5682-7520 


4558.90350-4740.54142 


120 


7660-9464 


4558.90350-4740.54142 
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Figure 1. Number of Ha spectra per bin of long and orbital phase. 



ing V2 r sin i = 60 ± 2 km s _1 . A new synthetic model for the donor 
was constructed using this rotational velocity. This model is used in 
the rest of this paper. The donor synthetic spectrum compared with 
a typical UVES spectrum reveals in general a good match except 
for minor deviations in short wavelength segments due to difficult 
normalization of the waving continuum (Fig. 2). 

Looking at luminosity-sensitive features like the Fe II & Ti 
II double blend at 4172-8 A, and similar blends at 4395-4400 A, 
4417 A and 4444 A, we do not observe evidence for systematic 
changes in luminosity for the donor star. The same is true for its 
temperature, as indicated by line strength ratios. 



3.3 Radial velocities for the donor 

Radial velocities for the donor were measured by cross-correlating 
the observed spectra with the FEROS spectrum taken at mid 
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Figure 2. The donor synthetic spectrum (DSS, dashed line) over-plotted with one of our UVES spectrum taken near main eclipse (<E> = 0.972, <J>; = 0.753). 
To get the DSS, the model synthetic spectrum has been scaled according to the light curve model (M12) to fit the percentage of light contributed by the donor 
at this epoch. It has also been shifted in wavelength to fit the donor radial velocity at this orbital phase. In the figure, the DSS is vertically shifted to match its 
continuum with the unity. Metallic lines were identified and labeled (A — 4500 A); uncertain identifications are marked with asterisks. 
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Figure 3. Radial velocities of the donor star (asterisks for CORALIE data 
and squares for UVES data) and the best fit for elliptical (up) and circular 
(down) orbits. Residuals are given in the inset graphs. 



Table 4. Donor radial velocities obtained by cross-correlation. Labels indi- 
cate UVES (U) and CORALIE (C) data. Typical RV error is 1 km s~ 1 . The 
complete table is given in electronic form. 



HJD 


RVfJons- 1 ) 


label 


2454558,90258 


-55,62 


U 


2454558,90352 


-54,56 


U 


2454558,90446 


-55,59 


u 


2454559,85331 


65,31 


u 


2454559,85425 


65,31 


u 


2454560,85597 


157,97 


u 


2454560,85691 


158,01 


u 


2454560,85784 


157,95 


u 


2454561,87559 


149,95 


u 


2454561,87653 


149,97 


u 



eclipse. After obtaining relative velocities with respect to this ref- 
erence spectrum, we applied a zero-point shift considering the tem- 
plate RV, measured by fitting the position of a set of known lines 
with simple gaussians and comparing these wavelengths with the 
corresponding laboratory wavelengths. The cross-correlation was 
performed in two regions deployed of H I and He I lines, viz. 4500- 
4800 A and 5050-5680 A. The radial velocities are given in Ta- 
ble 4. 

The RVs can be fitted with a sinusoid with radial-velocity half- 
amplitude K2 = 169.7 ± 0.6 km s^ 1 and zero point 7 = -19.6 ± 0.4 
km s ~~ 1 . These figures are significantly different from K2 = 1 8 1 ± 3 
and 7 = -2 ± 2 km s _1 obtained from the analysis of the single line 
Mgll 10917 by Ml 0, but this line is probably partly blended with 
Pa7 and He 1 10916 and affected by residual emission. A careful 
inspection of Fig. 3 shows non-random residuals for the circular fit. 

In order to resolve the question about the possible ellipticity 
of the orbit, we used the genetic algorithm PIKAIA developed by 
Charbonneau (1995) to find the orbital elements for V393Sco. 
The method consists in finding the set of orbital parameters 
that produces a series of theoretical velocities that minimize the 
function % 2 defined as: 

X 2 (P ,r,w,e,if2,7) = 



JV 

jv-6 / , y 
3=1 



V(tj,P ,T,ui,e,K2,i) 



(1) 



where N is the number of observations, P a is the orbital period, uj 
the periastron longitude, r the time of passage per the periastron, e 
the orbital eccentricity, K2 the half- amplitude of the radial velocity 
for the donor and 7 the velocity of the system center of mass. Vj 
and V are the observed and theoretical radial velocities at tj . The 
theoretical velocity is given by: 

V(t) = 7 + K 2 (cos(lj + v(t)) + e cos(u)) (2) 

where v is the angular velocity around the system center of mass 
obtained resolving the following two equations involving the 
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Table 5. Orbital elements for the donor of V 393 Scorpii obtained by mini- 
mization of the x 2 parameter given by Eq. (1) for the cases of circular (Case 
C) and elliptical (Case E) orbits . The value r* = r — 2454507.78 is given 
and also the error measured as half the difference between the maximum 
and minimum quantity in one isophote la. 



Case 


Po 

(days) 


T* 


e 




K 2 
(kms^ 1 ) 


7 

(kms" 1 ) 


x 2 


C 


7.7139 


9.6 


0.00 


1.95 


170.7 


-19.1 


9 


E 


7.7130 


1.3 


0.04 


1.49 


169.7 


-18.6 


2 


error 


0.0020 


0.1 


0.02 


0.08 


2.0 


2.8 





eccentric anomaly E: 

ton(|) = 
E-E sin(E) = §?(*-t) 



tan(§ ) 



(3) 
(4) 



A range of physically reasonable parameters need to be con- 
sidered so that the method works. For the period we used the range 
0-10 days, the eccentricity was set between and 1, uj between 
and 2tt, t between the minimum julian day and this value plus P a , 
K 2 between and (V max - V min ) and 7 between V min and V max . 

The most reliable way to get error estimates for this genetic al- 
gorithm is by Monte Carlo simulations, specifically by perturbing 
the best fit solution and computing the \ 2 of these perturbed solu- 
tions. To find the standard deviation region (a) encompassed by the 
joint variation of two parameters with all other parameters at their 
optimized values, we draw the contour corresponding to that value 
of Ax 2 for 2 de grees of freedom that includes 68.3% of the prob- 
ability. In our case this corresponds to A\ 2 = 2.30 (Bevington & 
Robinson 1992, Chapter 11, pp. 212). 

The best orbital elements for the circular and eccentric solu- 
tions along with their estimated errors are presented in Table 5. It is 
clear from the % 2 value in this table and the residual plots in Fig. 3 
that the elliptical solution provides the best fit, since it gives resid- 
uals without systematic trends and also the smaller \ 2 value. We 
note that our eccentric solution gives a small e value (0.042) but it 
is highly significant, according to the statistical test "pi" of Lucy 
(2005). 

It has been pointed out that gas stream and circumstellar mat- 
ter can distort spectroscopic features in semi-detached interacting 
binaries, producing skewed radial velocities and artificial small ec- 
centricities (e.g. Lucy 2005). For a non-interacting binary with the 
stellar and orbital parameters of V 393 Sco, dynamical tides should 
have circularized the orbit and synchronized the rotational periods 
(Zahn 1975, 1977). This should imply that the observed small ec- 
centricity is very likely spurious. However, the system is proba- 
bly found with a massive circumprimary disc remanent of a recent 
mass transfer burst and with the primary rotating at critical velocity 
(Ml 2). It is possible that the observed eccentricity could be the re- 
sult of a dynamical perturbation introduced by the circumprimary 
massive disc. This is an open question we let for a future investiga- 
tion. 



3.4 The He 1 5875 line and the mass ratio 

Radial velocities for He I 5875 were measured using simple 
gaussian fits. The corresponding spectral region is virtually uncon- 
taminated by donor spectral features. The RVs are in anti-phase 
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Figure 4. Up: He 1 5875 radial velocities at different long cycle phases, 0.2 
<$l^ 0.8 (circles) and 0.8 < s£ 0.2 (dots) and the best sinus fits. 



with the donor RVs, revealing their origin in or around the primary 
star (Fig. 4). There are large RV dispersion at a given phase, and 
the overall pattern of the RV curve changes during the long-cycle. 
This fact suggests that He I lines could be affected by residual 
emission at some epochs. It is therefore difficult to estimate the 
gainer radial velocity half-amplitude K\ and the mass ratio from 
these RVs. We notice, however, that the half-amplitudes are 82.9 
± 1.5 and 68.1 ± 4.5 kms _1 at low and high state, respectively 
(Table 6), much larger than K\ = 45.6 ± 4.6 km s _1 derived from 
the study of UV and IR lines by M10. If we assume that these 
half-amplitudes represent the motion of the more massive star then 
the larger values indicate a mass ratio q = M2/M1 « 0.4 or 0.5, 
whereas the M10 figure indicates q = 0.27 ± 0.03. M10 preferred 
the lower value for the mass ratio since: (i) it allows the formation 
of a disc around the gainer and (ii) it is more compatible with the 
mass derived for the secondary star. We add a third reason: (iii) a 
low mass ratio value is compatible with synchronous rotation for 
the secondary star. In fact, for a secondary star corotating with the 
binary: 



K 2 



■«) 



0.49g 2/3 
0.6 9 2 / 3 +ln(l+ 9 1 / 3 ) 



(5) 



(Eggleton 2006, Eq. 3.9), where K2 is the donor radial velocity 
half-amplitude. Using the above equation, K2 = 170 km s _1 and 
V2rSin i = 60 km s _1 , we obtain q w 0.27, consisting with the low 
q found by M10. 

The strength and shape of the He 1 5875 line vary during the 
orbital and long cycles. The FWHM increases systematically 
from $; = 0.5 to = 0.0 and decreases from high to low state. The 
FWHM increases during the orbital cycle attaining a maximum 
at $ = 0.9, similarly to He 1 10830 (M10). The amplitude of the 
FWHM variability is about 200 km s _1 . Emission at the flanks 
of the He 5875 absorption line is seen occasionally, with typical full 
peak separation AA5875 ~ 700 km s -1 . 

The long-cycle variability can be explained by filling emis- 
sion increasing toward the high state. At low state we observe the 
deeper and narrower line, and at high state the line is almost filled 
by emission, producing the observed FWHM increase. Accord- 
ingly, the equivalent width (EW) of this line anti-correlates with 
their FWHM. The maximum at $ = 0.9 is produced by an en- 
hanced red absorption wing. This event coincides with the passage 



A cyclic bipolar wind in the interacting binary V393 Scorpii 1 



0.80 
0.75 - 

0.70 
0.6a -| 
0.60 
0.55 - 
0.50 
0.45 - 
0.40 
0.35 - 
0.30 
0.25 



L3/(L1+L3+Ld) 




0.0 0.1 0.2 0.3 01 0.5 O.f 



0.9 1.0 1.1 1.2 1.3 



Figure 5. Light contribution factor for the donor at different spectral lines 
and orbital phases according to the M12 model. L\, L2 and are the 
gainer, donor and disc fluxes, respectively. 
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Figure 6. Equivalent widths and equivalent width ratios for synthetic spec- 
tra calculated for stellar atmospheres with logg = 4.0, V tur b= 2 km s _1 
and solar metallicity. 



of the hot spot across the line of sight. It could reflect large infall 
velocities in the absorbing material of the gas stream. 



3.5 Spectral disentangling and donor-subtracted spectra 

The variable and multicomponent nature of the spectra make diffi- 
cult their analysis and requires some kind of useful approximation 
to isolate single component contributions. In order to remove the 
donor light from the spectra we assume that its contribution to the 
total light is additive to the other light sources and represented by 
the model proposed by M12, i.e. it is independent of the long cycle. 
This assumption is justified because: (i) the amplitude of the long 
cycle is small compared with the large magnitude changes observed 
during the orbital cycle and (ii) we do not detect large systematic 
changes of the donor during the long cycle. 

We constructed donor-subtracted spectra by removing the 
synthetic spectrum of the donor from the observed spectrum. 
Firstly, the synthetic spectrum was Doppler corrected and scaled 
according to the contribution of the secondary star at a given 
orbital phase and at the given wavelength range. For normalized 
fluxes: 

f res (X, <S> , $ ; ) = /(A, $ , <J> ; ) - p(* , A c ) x f d (X, * ) (6) 

where f res is the donor-subtracted flux, / the observed flux, id 
the synthetic donor spectrum, p is the fractional contribution of the 
donor derived from our model and A c the representative wavelength 
where this factor is calculated. Theoretical p-factors correctly ac- 
count for the variable projected area of the donor and its larger flux 
contribution at longer wavelengths (Fig. 5). The method described 
above effectively removed the main contribution of the donor to the 
observed spectra. The result was a set of "donor- subtracted" spectra 
that were normalized to the new continuum. 

Disentangling of the gainer and the disc was not intended since 
there are no model available for the disc spectral distribution and 
probably its spectrum combine in a complex way with the spectrum 
of the primary. 



3.6 The He 1 4471 and Mg II 4482 lines 

The ratio between the equivalent width of He 1 447 1 and Mg II 4482 
is a good temperature indicator in early B-type stellar atmospheres, 
being larger for late B-type stars. We calculated equivalent widths 
and equivalent width ratios (R) as a function of stellar effective 
temperature for our grid of synthetic stellar spectra (Fig. 6). These 
calculations are useful to get some insight on the physical condi- 
tions of the region forming the helium absorption. 

EWaa71 and R show orbital and long-cycle variability. Dur- 
ing high state magnesium emission increases faster, producing a 
larger R at this state (Fig. 7). This result suggests the R could be 
used as a rough diagnostic for the long cycle stage in absence of 
a clear ephemeris, for old spectroscopic observations for instance. 
The orbital variability of R and EWavii at low state suggests that 
the effective temperature for the hot component varies with the or- 
bital phase and it is larger at 3> = 0.4. At low state EWmji and 
R are compatible with line formation in a photosphere with Thot 
= 16.000 K. This is the same temperature found for the inner disc 
region (and the gainer) by M12. In contrast, at high state R is larger 
than any theoretical value, mostly due to the weakness of the mag- 
nesium line by filling emission. This makes uncertain any i?-based 
temperature at this stage. 

3.7 Balmer and He 1 6678 lines 

3.7.1 General considerations 

According to the M12 model, at <E> =0.5 occurs the largest oc- 
cultation of the donor by the disc, that attains maximum visibility. 
We compare Ha and He 1 6678 A profiles at this phase for sev- 
eral long-cycle epochs and conclude that (Fig. 8): (i) at low state 
Ha shows a relatively narrow and deep central absorption (CA) 
flanked by emission peaks of similar intensity. CA is up to 65% the 
continuum level at this orbital phase, (ii) when the system goes to 
high state the global emission strength increases and the CA tends 
to disappear. After a careful examination of the whole dataset, we 
concluded that this behavior produces the variability of the peak 
separation described later, i.e., the stronger the CA, the larger the 
peak separation. Some flat-toped profiles observed at long maxi- 
mum fit this tendency. 
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Figure 7. He 4471 equivalent width and He4471/MgII4482 equivalent 
width ratio versus orbital phase. UVES data for 0.8 < < 0.2 (filled 
circles) and 0.2 ^ "3>; ^ 0.8 (open circles) are shown. 



We compare theoretical photospheric spectra for a critically 
rotating gainer with the observed Ha and He 1 6678 central absorp- 
tion profiles and conclude that they do not match at long cycle min- 
imum (Fig. 8). However, at high state the He I profile fits the theo- 
retical profile of a slower gainer. We propose that this result is only 
the consequence of the evident filling emission. In this view the 
broadening of 200 km s -1 probably reflects the kinematics of the 
gainer pseudo-photosphere. Others interesting features observed in 
these spectra are the Ha emission wings extending beyond 1000 
km s _1 (see Section 4.2). 

The general variability for these lines is illustrated in Fig. 9. 
The helium line is filled by emission at high state, a characteristic 
also seen in other helium lines. It is notable the presence of emis- 
sion flanks with AAe678 ~ 900 km s~ . These flanks are also seen 
with varying intensity in other helium lines. 

The Ha line appears as a deep absorption surrounded by weak 
emission shoulders during low state but it is a strong and double 
emission (or in few occasions triple or irregular) during high state. 
For double emission lines the violet emission peak changes its rel- 
ative intensity with respect to the red peak during the orbital cycle. 
This property is measured through the V/R ratio, where V is the 
violet peak intensity and R the red peak intensity. The V > R 
condition near $ D = 0.75 reveals a larger Ha emitting region at 
the binary hemisphere containing the hot spot (Fig. 9). Note also 
in Fig. 9 the presence of a weak emission line redward Ha that we 
identified with C 1 6588. This line follows the donor RV curve and 
will be discussed in Section 3.9.2. 

The H/3 and H7 emission profiles show similar behavior to 
Ha, but with weaker emission and most importantly with very ex- 
tended emission wings (EEWs), beyond ± 2000 km s _1 (Fig. 10). 
These velocities are too large to be explained in terms of Keple- 
rian rotation around the gainer. The EEWs usually appear at the red 
profile flank during the first half of the orbital cycle and at the blue 
flank during the 2nd half, consisting with an emitting region in the 
binary hemisphere containing the hot spot, a region high enough to 
be observed as an emitting source rather than an absorbing media. 
As occurs with the low-velocity B aimer line emission, the strength 
of these EEWs is larger at high state. 

Occasionally we observe a narrow emission feature redward 
He 1 4921 that is sometimes accompanied with a redshifted absorp- 



1.2 
1.1 
1 

0.9 
0.8 
0.7 
0.6 



45. 07. 3 
0.55, 25, 2 



^\^ a49 " 5 4B 0^5.'3 
iJTMV n 5n n 77 3 



50. 77. 3 
vsini = 415 km/s 




-500 500 

velocity (km/s) 



1,3 



■ 0.917(5)+0.279(8) sin[n(x+0.1 8(1 ))/0.502(4)] | 



1,0 
0,9 



s 
x 



0,6 



i 


i 1 i ■ i 1 i ■ i 

15 

1 2 ,. 








/ * \ o 5 








/ \ ° 




1 " 






0,2 0,4 0,6 0,8 

orbital phase ■ 




i 







0,0 0,2 0,4 0,6 0,8 1,0 1,2 1,4 1,6 1,8 2,0 

long-cycle phase 





0.45. 0.07. 3 
55. 25 2 
49-0 54.0.34. 11 








0.48, 0.65, 3 
50, 77, 3 
vsini = 200 km/s 
























III 




























HP 










































m I.JI 











































































































-500 500 

velocity (km/s) 



Figure 8. Up and down: Donor-subtracted Ha and He I 6678 profiles near 
secondary eclipse at different long cycle phases. Labels indicate, from left 
to right, orbital phase, long cycle phase and number of spectra averaged. 
Synthetic profiles of the gainer calculated with parameters given by Ml 2 
are shown for different rotational velocities. X-axis velocities are with re- 
spect to the system center of mass. Middle: The behavior of the intensity 
of the central absorption of the donor-subtracted Ha profile at secondary 
eclipse during the long cycle. The best sinus fit is also shown along with the 
histogram of orbital phases for the considered spectra. 
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Table 6. Results of the sinusoidal fits (7 + K sin(2-7r(a; — 5))) to the central 
absorption RV curves. The root mean square of the fit is also given. The 
parameters 7, K and rms are given in km s~ 1 . 



Line (state) 


7 


K 


6 


rms 


Ha (all) 


-7.9 ± 1.5 


59.6 ± 2.3 


0.550 ± 0.005 


20.1 


H/3 (low) 


-9.8 ± 1.0 


64.0 ± 1.5 


0.546 ± 0.004 


15.1 


H/3 (high) 


-4.1 ± 3.5 


78.4 ± 5.4 


0.525 ± 0.009 


32.1 


H7 (low) 


-15.8 ± 1.1 


71.3 ± 1.5 


0.529 ± 0.003 


14.7 


H 7 (high) 


-14.8 ± 2.7 


83.2 ±4.4 


0.517 ±0.006 


26.6 


OI 7773 (low) 


-13.5 ±2.6 


32.1 ± 3.8 


0.593 ± 0.017 


20.7 


He 1 5875 (low) 


-12.3 ± 1.1 


82.9 ± 1.6 


0.572 ± 0.003 


15.3 


He 1 5875 (high) 


3.0 ±3.1 


68.1 ±4.5 


0.555 ± 0.010 


25.0 


He 1 6678 (all) 


-8.8 ± 1.3 


70.2 ± 1.9 


0.554 ± 0.004 


21.6 



tion (Fig. 10). We identify this feature as CI 4932; its RV suggests 
an origin in or around the secondary star. This line will be discussed 
along with other weak spectroscopic features in a forthcoming pa- 
per. 

3.7.2 Radial velocities of emission line features 

Here we present the analysis of the radial velocities of the cen- 
tral cores and emission peaks of the Balmer lines and some he- 
lium lines. The CA RV curves were fitted with sinus functions with 
the best fitting parameters given in Table 6. In general, the radial 
velocities for Balmer and helium features are in anti-phase with 
the donor RVs, indicating an origin somewhere near the gainer 
(Figs. 11 and 12). However, they do not follow the expected mo- 
tion of the gainer. The half-amplitudes are usually larger than the 
adopted K\ value and the positions of the light-centers, as mea- 
sured by the S parameters, are displaced from the center of the pri- 
mary. We notice larger Balmer amplitudes during high state and the 
opposite for Helium lines. We also notice the relatively small am- 
plitude of the 1 7773 line. The line centers, as measured by the 
7 parameter, have larger velocity than the systemic velocity, which 
can be interpreted in terms of an outflow or wind. 

Contrary to He 5875, the He 1 6678 CA shows small changes 
through the long cycle. The RV of the blue emission peak shows 
a quasi-sinusoidal behavior with large amplitude, whereas the 
red peak remains almost stationary at 390 km s _1 shifting to 550 
km s _1 at high state. On low state the blue Ha emission peak shows 
a large blueshift at $ =0.25, as occurs with He 1 5875 (Fig. 4). 
This might be due to the approach of the bright spot at longitude 
162° (M12). Another notable feature is the shift to bluer wave- 
lengths of the Ha CA around $ = 0.15, as corroborated by data 
obtained on high state at 3 different observatories. This anomaly 
is also observed in the asymmetric absorption wings discussed in 
Section 4.5. 

If we use the system inclination i — 80° (Ml 2), and the basic 
kinematics formula: 



where r represents the radial distance from the center of rotation, 
P the rotational period in days and v the linear orbital velocity of 
material moving in the orbital plane measured in km s _1 , then we 
can find the position of the light-centers of the lines listed in Ta- 
ble 6. We find all line light centers located practically in the line 
joining the stellar centers. In addition, light centers for Balmer and 
helium lines are located between 9 and 13 Rq from the barycen- 



ter measured in the direction of the gainer; for 1 7773 this figure 
is 4.9 Rq. Remembering that the center of mass is located at 7.16 
Rq from the gainer and that this star has a radius of 4.4 Rq (M12), 
these results show that all lines are formed somewhere around the 
gainer, 017773 probably near the gainer hemisphere facing the 
donor (consistent with the location of the DACs forming region) 
and the other lines in more distant regions. Since the gas dynamics 
is very likely not limited to the orbital plane and includes impor- 
tant vertical motions not considered in this analysis we should take 
these results with caution. 



3. 7. 3 Equivalent widths and peak separations 

On low state the equivalent width of the He 1 6678 line shows large 
orbital variability and increasing emission during the main eclipse 
(see Fig. 11, similar behavior is observed in HeI4471 in Fig. 7). 
This simple observation reveals that the source of the line emis- 
sion is not eclipsed as occurs with the continuum flux of the donor- 
subtracted spectra. Interestingly, the same parameter shows much 
weaker orbital variability on high state, and importantly the strong 
variability during eclipse disappears. This behavior is difficult to 
understand, but it could be partly explained if the extra continuum 
light present at the high state is not eclipsed. In other words, the 
wind not only contributes to line emission at low/high state but also 
to continuum emission at high state, i.e. it is optically thicker during 
the long maximum. 

As mentioned before, the Ha peak separation is smaller at 
high state, principally due to the weaker CA. Moreover, the equiv- 
alent width of Balmer lines vary during the long cycle in anti- 
correlation with the peak separation (Fig. 13). The Ha peak sepa- 
ration goes from 690 at low state to 340 km s _1 at high state, where 
it cannot be due to Keplerian motions since it should imply an orbit 
outside the gainer Roche lobe. For H/3 and H7 the corresponding 
figures are 800 to 500 km s _1 and 1000 to 650 km s~ , respectively. 
This behavior is opposite to the observed in He 6678, where the 
peak separation increases on high state (Fig. 10). 

If the CA is formed in the gainer disc-like extended photo- 
sphere, then we can reasonably assume that rotational broadening 
determines the CA width, and also indirectly the peak separation. 
The maximum rotational velocity in the disc is derived from AA 7 
at low state (about 500 km s _1 ) and suggests that the primary is 
in fact rotating at critical velocity. As for helium lines, the Balmer 
emission increases during main eclipse, indicating that the main 
source for Balmer line emission is not eclipsed neither. We mea- 
sure maximum equivalent widths at main eclipse of -40, -20 and 
-20 A, for Ha, H/3 and H7, respectively (Fig. 14). 

3.7.4 Ha V/R variability 

The Ha Vj R ratio follows a quasi-sinusoidal pattern during the or- 
bital cycle with local maxima at $ D = 0.05, 0.45 and 0.70 (Fig. 15). 
In the same figure a reference sinusoid represents the smoothed 
orbital V/R variability. This variability can be explained by the 
displacement of the CA across a stationary emission profile. This 
view is difficult to test measuring the supposedly stationary pro- 
file extension, due to their strength variability and gradual merging 
with the continuum, but it is supported by the fact that the V < R 
branch coincides with the blue-shifted CA position and the V > R 
one with the red-shifted. The stability of the V/R curve at high and 
low state reflects the stability of the disc where the CA is formed, a 
result also supported by our light-curve photometric model (M12). 
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Figure 9. Donor-subtracted Ho (left) and He 1 6678 (right) spectra at se- 
lected epochs. Orbital phases (left) and long-cycle phases (right) are given. 
Note the C 1 6588 emission redward Ho. 



Figure 11. Up: Radial velocity for the central absorption and the red and 
blue peaks of the He 1 6678 line. The best sinus fit for the central absorption 
is also shown. Down: The equivalent width for the He 1 6678 as a function 
of the orbital and long cycle phase. The color band represents <E>; ranges. 
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Figure 10. Donor subtracted spectra showing H/3 extended emission wings 
up to ± 2500 km s -1 . The He 1 4921 absorption line is also visible along 
with narrow C 14932 emission. The numbers between parenthesis indicate 
the number of averaged spectra. 

The zero point of the V/R curve is larger than unity, indicating 
the predominance of the V > R condition. This is consistent with 
an outflow of the emitting region, as also indicated by the 7 values 
in Table 6. The V/R curve shows 3 enigmatic successive and rapid 
reversals just after main eclipse. In addition, a "s-type" variation is 
found around the referential sinusoid between <E> = 0.7-1.2 and a 
"z-type" perturbation between <E> = 0.4-0.7. 

In order to interpret the s-type and z-type waves, we assume 
that the high latitude region responsible for Balmer emission ex- 
tends to the disc of such way that the lower parts merge with the 
disc bringing its rotational velocity. In this case the s-wave can be 
viewed as the eclipse of the approaching disc edge by the donor 
during $ =0.7-1.0 and the eclipse of the receding edge during 
<J> = 0.0-0.2. In the same picture the z-wave probably shows the 
eclipse by the disc of the receding (approaching) hot spot region 
in <I> = 0.4-0.5 (0.5-0.6). We remember that according to the light 
curve model a hot spot is expected in the place where the gas stream 
hits the disc. This hot spot faces the observer at <!> = 0.9 (M12) and 



probably produces, as already mentioned, the FW H M 5 ^^ peak at 
$ o =0.9. 

In summary, our observations seems to indicate that a wind is 
emerging from the disc producing Balmer emission lines by elec- 
tronic recombination. We notice that while our wind-disc model 
explains well the Balmer line variability, the alternative hypothesis 
of a circumbinary disc cannot explain many of the observed fea- 
tures. 

3.8 Balmer difference spectra 

In order to analyze the variable line emission we constructed dif- 
ference spectra by subtracting profiles from a reference profile ob- 
tained roughly at the same orbital phase but different long cycle 
phase, usually that one with the minor line emissivity. We do not 
correct for the underlying variability of the continuum through the 
long cycle, since it should have a minor effect on line profile in- 
tensity. Difference spectra constitute useful representations of line 
emission changes at a given orbital phase and projected velocity 
bin during the long cycle. The results for Ha discussed here are 
basically the same that for H/3 and H7. 

All difference spectra for Ha taken around $ = 0.0 indicate 
that the emission grows principally at low projected velocities when 
the system approaches the high state (Fig. 16). The characteristic 
single peak profile of the extra emission does not correspond with 
disc emission, either circumprimary or circumbinary and reinforces 
the thesis of a wind. 

At $ = 0.75 the emission grows around -200 km s _1 and 50 
km s _1 (Fig. 16). These velocities suggest that emitting material 
exist around the donor and the primary, a fact corroborated by the 
Ha Doppler map later. In general, Balmer difference spectra show 
that the wind line emissivity increases on the high state. 

3.9 Doppler Tomography 

Doppler tomography was introduced as a tool for the study of ac- 
cretion discs by Marsh & Home (1988). This analysis method is 
now a widespread procedure in the study of emission lines in Al- 
gols (Richards 2004), and provides a quantitative mapping of opti- 
cally thin line forming regions in the velocity space. When severe 
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Figure 12. Up: Radial velocity for the blue and red peak of the Ho emis- 
sion and the central absorption. The best sinus fit (excluding outliers around 
$o =0.15) along with the assumed track for the gainer are also shown. 
KPNO data are shown as asterisks and are not included in the fit. Middle 
and Down: The radial velocity for the H/3 (middle) and H7 (down) central 
absorption. The best sinus fits on low state are shown, along with the as- 
sumed track for the gainer. In all graph data for 0.8 < < 0.2 (max) and 
0.2 <&; ^ 0.8 (min) are shown with separated symbols. 



Figure 13. The peak separation (up) and equivalent width (down) of the 
donor-subtracted Ha line profile. Data for 0.95 < <E> < 1 .05 are not shown. 
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Figure 14. The equivalent width of the donor-subtracted Balmer line pro- 
files during main eclipse. Data on the long cycle phases 0.8 < <&; < 0.2 
are shown as filled circles and those with 0.2 ^ "3?; 0.8 as open circles. 
Asterisks indicate KPNO data. 



self-absorption and/or intrinsic line broadening are present then the 
tomograms provide at least a concise and convenient way of dis- 
playing phase-resolved line profile measurements. The Doppler re- 
constructions presented in this paper were computed using the fil- 
tered back-projection method (Rosenfeld & Kak 1982). The two- 
dimensional reconstruction of structures in the orbital plane has a 
unique solution that depends only on the two-dimensional data set 
(i.e. f v ($ , A) where f v is the flux). Our code has been success- 
fully applied in a number of cases (e.g. Mennickent, Diaz & Arenas 
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subtracted Ha double emission profiles (the inset shows a zoom around the 
main eclipse). The sinusoid V/R = 1.03 - 0.15 sin(27r(<E> o -0.05)) has been 
plotted as reference. Data with 0.8 < <J>; < 0.2 (max) are shown as filled 
circles and those with 0.2 ^ <I>; ^ 0.8 (min) with open circles. Asterisks 
indicate KPNO data. 
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1999) and the comparison of observed and reconstructed line pro- 
files shows the self-consistency of the method. 

When interpreting Doppler maps we should keep in mind 
the difficulty in the representation of optically thick structures and 
components above or below the orbital plane. This is specially rel- 
evant for the Balmer emission lines discussed in this paper. The 
evidence presented in this section for Ha is not conclusive, and the 
map should be only taken as a particular two-dimensional repre- 
sentation of the three-dimensional mass flows in the binary. 



3. 9. 1 The Ha Doppler map 

We built the Doppler map for the donor-subtracted Ha line with 
resolution FWHM = 65 km s _1 . Only spectra with 0.8 <<!>;< 
0.3 were used, since others showed too weak emission. The map 
is shown in Fig. 17. The stream path plotted in the map is for a 
point mass primary. For a 4.4 Rq primary radius the stream would 
ultimately collide with the star at tomogram velocities V x = -675 
km s _1 and V y = -310 km s _1 (i.e. in the lower left quadrant) or 
even sooner with disc or wind. 

The Ha absorption profile goes below the continuum during 
some phases. This is inconsistent with an optically thin emission 
formed in the gas velocity field. Having said that 3 interpretations 
for the Doppler map remain: (i) the gas is optically thick in the 
lines, (ii) the photospheric spectrum has a broad absorption com- 
ponent and (iii) a combination of both is present (an optically thick 
wind and/or disc). The map indicates that the absorption has a 
FWZI of circa 150 km s _1 and has significant radial velocity 
with phase still consistent with the primary. We computed phased 
spectrograms and the orbital velocity modulation of the absorption 
can be confirmed at high state. It is possible to interpret this region 
as the optically thick circumprimary disc. Additional pieces of in- 
formation are that the emission almost disappears at L\ and around 
$o = 0.65. Finally, emission is detected around the donor, confirm- 
ing the results inferred from our Ha difference spectra. In sum- 
mary, and taking into account our complementary spectroscopic 
analysis, the Ha Doppler map is consistent (but not conclusive) 
with an optically thick circumprimary disc plus high-latitude wind 
and emitting material located around the donor star. 
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Figure 16. Up: Ha difference spectra near <J> = 0.00 at selected long cycle 
phases. Down: Ha difference spectra near <I> = 0.75 at selected long cycle 
phases. Velocities are with respect to the system center of mass. 



3. 9. 2 The Mg II 4481 and C 1 6588 Doppler maps 

We detected Mg II 4481 and CI6588 emissions in the donor sub- 
tracted spectra whose velocities indicated an origin around the sec- 
ondary star. These emission are very variable in intensity and some- 
times show double peaks. They are present at high and low state. A 
renormalization of the continuum was made around these lines and 
Doppler maps were constructed in the usual way. 

The Doppler map for Mg II 4482 was centered at 4481.2 A, it 
shows that this line is formed in the donor and it is not produced by 
irradiation of the hemisphere facing the gainer, but probably cor- 
responds to chromospheric activity (Fig. 18). The same result is 
obtained with the CI 6588 line. The presence of near/underlying 
absorptions in this line and the data sampling and spectral noise 
produces strong streaks mimicking emission at L2, L3 and the gas 
stream region. These are very likely artifacts of the reduction pro- 
cess. 

The discovery of chromospheric emissions open the possibil- 
ity that the donor in V 393 Sco posses a non-neglected magnetic 
field. 
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Figure 17. The back-projection Doppler map for the donor- subtracted Ha 
line during high state (0.8 < <!>; < 0.3). Crosses show, from up to down, 
donor center of mass, LI point, system center of mass and gainer center of 
gravity. Roche lobe surfaces for a point mass gainer are shown along with 
the gas stream path (lower track). The upper track represents the keplerian 
velocity of a disc if it would exist along the stream. 



3.10 The 1 7773 triplet and the Si II 6347-6371 doublet 

The 017773 line appears in absorption on low state and al- 
most completely filled by emission on high state. This line is 
a triplet with the redder component displaced by 67 km s _1 and 
107 km s _1 from the two other components, but in our spectra it 
appears broadened and the single components are not resolved. 
The line shows pronounced blue absorption wings in the phase 
range <1> = 0.55-0.84 and enhanced red absorption wings during 
$ = 0.20-0.45 (Fig. 19). These line depressions are similar to 
those observed in He 1 10830 (M10) and are present even before 
applying the disentangling process. 

We detected Discrete Absorption Components (DACs, of typ- 
ical depth 3 to 5 % of the continuum intensity) in virtually all 
017773 spectra with enhanced absorption wings, but in much 
larger number in the blue wing. DACs are visible even before spec- 
tral disentangling and in spectra free of fringes and bad pixel re- 
gions. From the above, and from their visibility at specific orbital 
phases, we argue that DACs are not artifacts of the disentangling 
process or instrumental setup. 

We averaged consecutive spectra taken few minute apart on a 
same night in order to improve the signal to noise ratio (typically 
S/N ~ 120). Then we measured DAC radial velocities with re- 
spect to the reference wavelength A = 7773.83 A. Up to 8 DACs 
were detected in a single spectrum but more could be present hid- 
den by the waiving continuum. Most DACs are observed with pos- 
itive velocities during the first half of the orbital cycle but with 
negative velocities the second half (Fig. 20). We notice that DACs 
velocities are scattered around the line representing the donor star 
velocity, and are profusely observed at blue and red wings around 
$o = 0.4. This is illustrated in the asymmetric distribution observed 
in DAC radial velocities (Fig. 20). The corresponding distribution 
of velocities relative to the donor is smooth and relatively symmet- 
ric. 

Interestingly, we noticed that there is no indication of the 




Figure 18. Up: The back-projection Doppler map for the donor-subtracted 
Mg II 4482 line. Symbols are as in previous figure. The FWHM resolution 
is 60 km s -1 . Down: Same as above but for the CI 6588 line with resolution 
40kms _1 . 



triplet structure in any of the DACs. In addition, 20% of the DACs 
are double, with separation between the components of about 45 
km s _1 . DACs have velocities a few hundreds of km s _1 relative 
to the binary center of mass and are detected at low and high state. 
It is notable that profiles taken 69 days apart, at <J> = 0.68 and 
= 0.50 (HJD 2454672.52566 and 2454703.49230) show roughly 
the same set of DACs at very similar velocities. At high state dis- 
crete emission features appear in the spectra and the identification 
of DACs is difficult. One is tempted to say that some DACs are 
replaced by discrete emission features at the high state. 

Stationary or almost-stationary DACs are observed in the red 
wing of the line, with velocities 265 and 365 km s _1 . These fea- 
tures could correspond to Fell 7780.4 and Fell 7783.0 with ve- 
locities 12 km s~ x relative to the bary center, but we cannot dis- 
card OI "satellite" lines (Fig. 20). Around $ = 0.2-0.3 the DAC 
at 365 km s _1 (eventually Fell 7783.0) attains the larger velocity 
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Figure 19. Up: Examples of donor-subtracted OI 7773 lines showing DACs 
in enhanced blue wing. Middle: Examples of donor-subtracted OI7773 
lines showing DACs in enhanced red wing. Down: Donor-subtracted 
01 7773 spectra during the high state. Labels indicate orbital and long-cycle 
phases. 



(470 km s _1 relative to the 1 7773 center) and is wider than in 
other epochs. In contrast to these almost-stationary features, we 
notice that some DACs seems to follow, at least at some epochs, 
quasi-sinusoidal patterns following the donor orbital motion. 

The RV of the 1 7773 CA presents a slight orbital oscilla- 
tion with half-amplitude 32 km s _1 in anti-phase with the donor 
and center displaced by 5 km s _1 from the systemic velocity. This 



behavior changes drastically at long cycle maximum, as seen in 
Fig. 20, adopting a saw-teeth pattern, with a large velocity (+125 
km s _1 ) prior to $ D = 0.5 and a low velocity (-125 km s -1 ) just af- 
ter. Difference profiles shown in Fig. 21 indicate that this behavior 
is the consequence of the disappearance of the red (blue) emission 
wing during <E> = 0.45 (0.59). This phenomenon can be interpreted 
as the occultation of wind high- velocity regions by the gainer or the 
disc. 

We also find DACs in Sill 6347/6371. The 50% of the DACs 
detected in Si II 6371 are also detected in Si II 6347. In addition, 
44% of the Sill DACs have RV correspondence with those ob- 
served concurrently in 1 7773. We find relatively strong station- 
ary lines at 6375.8 A and 6379.2 A that we identified with diffuse 
interstellar bands (Herbig 1975). 

Radial velocities of broad donor-subtracted Sill 6347/6371 
line profiles present large scatter, but they are roughly in anti-phase 
with the donor RVs revealing their origin in or around the gainer. 

3.11 Enhanced absorption wings 

Enhanced absorption wings are observed in several lines, especially 
1 7773 (Fig. 19, upper panel) and He 1 5875 (not shown). We mea- 
sured the velocity of the wing edge with respect to the line center 
resulting in averages of 350 ± 25 (standard deviation) and -350 ± 
20 km s _1 for the red and blue wing of 1 7773 and 410 ± 55 and 
-390 ± 50 km s -1 for the red and blue wing of He 1 5875, respec- 
tively (Fig. 22). We propose that enhanced absorption wings are 
caused by photon absorption in the gas stream and circumprimary 
disc. 

Blue depressed wings occur in the range <1> = 0.5-0.9, roughly 
coinciding with the phases when the stream faces the observer, 
and red enhanced wings in the range <3> = 0.9-0.45, just when 
the observer see the receding stream (the exceptions are the blue 
wings around $ = 0.15 and the OI red enhanced wings already 
at <I> = 0.85). The blue wing at $ D = 0.15 could be explained if 
part of the stream bounces back after impacting the disc produc- 
ing blue-shifted additional absorption. This behavior is similar but 
not strictly the same that the one followed by DACs. In particular, 
the red enhanced wing appears around <J> = 0.9 but at this phase 
DACs still are present in the blue wing. Another interesting result 
is the abrupt change between positive and negative velocity, with 
no smooth transition between them. 



4 DISCUSSION 

Many of the phenomena described in this paper can be understood 
in terms of a global mass-loss process involving the transfer of 
mass between the donor and the gainer. In the following we discuss 
the spectroscopic variability in terms of possible line forming re- 
gions, line forming mechanisms and mass loss/exchange channels. 
A comprehensive summary of all these issues is given in Table 7. 



4.1 The wind revealed in the helium emission 

It is unusual to find in a binary star stationary spectral features, 
since most of the emitting regions should bring the fingerprints of 
the orbital motion. Therefore the discovery of stationary red emis- 
sion in He 1 6678 is important to constrain the physical scenario. 
The first thing we notice is that the RV variability of the emission 
peaks cannot be attributed to the motion of the central absorption 
feature, as in the case of Ha. This feature is in the middle of both 
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Figure 20. Radial velocities refereed to the binary center of mass for DACs (circles) and line center (filled circles for 0.39 <<!>;< 0.76, filled squares for 0.05 
<<!>;< 0.07) for the OI7773 line in UVES spectra. Measurements for Sill 6347 DACs (triangles) and Sill 6371 DACs (inverted triangles) are also shown. 
The theoretical radial velocity for the donor is shown by a sinusoid and dashed lines of constant velocity 0, 265 and 365 km s — 1 are also displayed. 



Table 7. Summary of spectroscopic variability, possible line forming regions and line formation mechanisms. ABS, CA, EM, RV and SIN stand for absorption, 
central absorption, emission, radial velocity and sinusoidal, respectively. Bound-bound emission is labeled b-b. 



Feature 


Variability 


Main Region 


Mechanism 


H I, He I, Si II & O I low velocity EM 


Larger at high-state 


wind 


b-b EM 


HIandHeICA depth 


Larger at low-state 


disc 


Self-ABS 


HI, He I, Sill&OICARV 


Quasi-SIN, low amplitude 


disc 


Self-ABS 


Balmer V/R curve 


SIN + perturbations on eclipses 


wind + disc 


Doppler + Eclipses 


He I V&RpeakRVs 


V quasi-SIN, R stationary and shifts at high state 


wind 


b-b EM 


DACs 


Mostly in depressed blue wings on 2nd half-cycle 


stream? 


unknown 


Balmer extended EM wings 


Red side 1 st half-cycle, blue side 2nd half-cycle 


wind 


photon scattering 


Enhanced ABS wings 


Complex, mostly blue side 2nd half-cycle 


stream + disc 


ABS 



peaks and seems not influence the RV of the stationary red peak. 
The more simple explanation for this behavior is assuming an al- 
most stationary broad single emission peak divided by the CA, and 
that this CA has an important blue wing component perturbing the 
measurement of the blue emission peak, especially at $ = 0.40, 
when the perturbation is large. The blue depressed wing could be 
due to the gas stream facing the observer at this phase and be the 



low latitude manifestation of the wind producing the broad helium 
emission. This wind could be produced by stream deflection at 
the disc/stream interaction region. Following the same reasoning 
that for Balmer lines, the helium half-peak separation (up to 450 
km s _1 ) is also consistent with a critically rotating gainer unable to 
accrete more material. 

In principle the helium emission wings could be due to scat- 
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Figure 21. Donor-subtracted OI7773 difference spectra during epochs of 
maximum RV shift at long maximum (<I> = 0.45 and 0.59) and zero RV 
shift ($ = 0.20). Labels indicate orbital phase and the difference between 
the long-cycle phases. 
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Figure 22. The velocity for the edge of the enhanced wing for 1 7773 and 
He 1 5875 relative to the line center. Averages values for 1 7773 are shown 
by dotted lines. Color bands indicate <3>; ranges. 



tering and true absorption both in the line and continuum by the 
Schuster's mechanism (e.g. Gebbie & Thomas 1968, Israelian & 
Nikoghossian, 1996). However, it is difficult to explain the asym- 
metric behavior between blue/red peak under this hypothesis. 

The presence of low velocity emission in He I lines confirms 
our previous view that the wind also emits in helium. As happens 
for Balmer lines, larger helium emission during the high state im- 
plies an increase in wind emissivity during this epoch. 

van Rensbergen et al. (2008) claim that gainer's spin-up and 
hot spots can drive mass out of a binary. They argue that releas- 
ing of accretion energy produces enough energy to sustain a wind 
from a small stream impact region. Their model gives a physical 
basis for a mass loss larger than the one driven by winds of isolated 
single stellar components of V 393 Sco. In addition, the discovery 
of a chromospherically active donor suggests that magnetism could 
play a role as a wind driver mechanism. The symmetry of the binary 
configuration regarding the orbital plane indicates that this wind 
should be bipolar. 



4.2 Balmer extended emission wings and the stream/disc 
interaction 

The presence of extended emission wings (EEWs) in some Balmer 
lines beyond ± 2000 kms _1 cannot be attributable to emission 
from a Keplerian disc. We notice that the wind projected veloci- 
ties (hundreds of km s -1 ) are too low to explain EEWs by elec- 
tronic recombination, and that the mechanism of Thompson scat- 
tering could produce extended emission wings as occurs in Be stars 
(Poeckert & Marlborough 1979) and the luminous blue variable 
AG Car (Schulte-Ladbeck et al. 1994). EEWs appear in the red 
wing on the 1st half-cycle and in the blue wing on the 2nd half- 
cycle so the forming medium shares the same orbital variability that 
the medium responsible for DACs and enhanced absorption wings. 
However, we notice that the emission features require a medium at 
latitudes higher than absorbing material. This suggests that part of 
the stream is deflected by the disc towards high latitudes produc- 
ing high density regions responsible for the Thompson scattering 
and EEWs. All of this activity occurs principally at the 1st and 2nd 
binary quadrant. 

4.3 Optically thick disc, hot spot and gas stream 

The central absorption cores of the lines studied in this paper are 
probably produced in the gainer pseudo-photosphere, modeled by 
Ml 2 as an optically- thick disc. The temperature inferred for the 
disc from spectroscopic line ratios T^ot = 16.000 K coincides with 
those inferred from the light curve analysis. The orbital behavior of 
line cores is consistent with light-centers placed somewhere along 
the line joining the stellar centers near the gainer. Variations in 
the depth of the CAs are explained by changing amounts of fill- 
ing emission as the wind increases its strength at high state. 

The gas stream is revealed at <& = 0.9 in the He I FWHM 
maximum and in the asymmetries observed in O I and He I lines 
(Section 3.11). The hot spot, on the other hand, shows up in the 
"z-wave" of the Ha V/ R variability. 

4.4 Line emission: the optically thin high latitude wind and 
donor activity 

The increase of line emissivity during eclipses and the sustained 
emission increase at low velocities during high state strongly sug- 
gest that most of the low velocity Balmer emission comes from 
electronic recombination at the high latitude wind. The same occurs 
for the emission seen in 1 7773 and the Si II doublet. As discussed 
in Section 3.7.2 part of the Ha emission also comes from the disc 
outer edge and a region around the donor. It is remarkable that when 
the system approaches the high state the line emission increases in 
the wind (the low velocity region) and around the donor (Fig. 16). 
This finding suggests a connection between the wind strength and 
the activity responsible for the line emission around the donor. This 
possible connection will be explored in a forthcoming paper dedi- 
cated to the study of weak emission features detected in the resid- 
ual spectra. So far, Doppler maps show that the Mg II 4481 and 
CI 6588 residual emissions are formed in a chromospherically ac- 
tive star. This finding arises the possibility of a magnetically active 
secondary in V 393 Sco. 

4.5 The DAC/stream connection 

We find that DACs follow the pattern of the enhanced absorption 
wings, so a connection between DACs and the gas stream is possi- 
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ble (see Sects. 3.10 and 3.11). The narrowness of DACs (FWHM 
a few tens of km s _1 ) indicates that they are formed in a differ- 
ent medium that the main 1 7773 line. Here we propose, tenta- 
tively, that DACs can be identified with tracers of donor mass loss. 
Traditionally, a continuum mass stream is considered the output of 
Roche-lobe overflow. However, the presence of DACs could reveal 
a stream composite of discrete and relatively narrow mass flows 
following magnetic field lines connecting the donor with the disc. 
In this view, DACs could be due to Zeeman splitting of absorption 
lines produced in streams projected against a bright continuum. The 
presence of discrete emissions on high state could indicate a re- 
ordering of the streams, carrying material by high latitudes produc- 
ing emission rather than absorption lines. The orbital distribution of 
DAC RVs could be result of the complex structure of the magnetic 
field. 

In order to support the above interpretation we notice that: (i) 
Zeeman splitting has been reported in the Si II doublet of the re- 
lated interacting binary j5 Lyrae (Skulskij & Plachinda 1993), (ii) 
the transient radio nature of V 393 Scorpii has been interpreted in 
terms of a stellar coronal magnetic field B ~ 10 G (Stewart et 
al. 1989), (iii) our discovery of a chromospherically active donor 
suggests a magnetic donor star and (iv) the discovery of a large 
and stable coronal loop observed in Algol and RY Vul (Peterson et 
al. 2010, Richards et al. 2010) indicates that evidence for magnet- 
ically channeled streams has been found in DPV-related interact- 
ing binaries. From the theoretical point of view, the physical basis 
for the donor magnetism might be its rapid synchronous rotation, 
that might generate a dynamo-powered magnetic field. This idea 
has been explored and vindicated for Algols with late-type com- 
ponents (e.g. Sarna et al. 1998, Moss & Tuominen 1997) and for 
rapidly rotating A-type stars (Featherstone et al. 2009). 

At present the nature of DACs is obscure, but our findings 
provide strong observational constraints for any competent model, 
and suggest that spectropolarimetry could be a successful technique 
in the diagnostic of the feature's splitting. 

4.6 New insights on the long photometric cycle 

In the past sections we have presented evidence for the existence of 
several circumstellar regions in V 393 Sco; an optically thick disc- 
like pseudo-photosphere, a gas stream, a hot spot, a chromospher- 
ically active secondary star and a high-latitude wind revealed in 
almost stationary Balmer and He I emission. We propose that all 
these regions and their associated phenomena indicate the presence 
of a mass-loss phenomenon in V 393 Sco that is cyclically mod- 
ulated with the long period. In this section we explore possible 
causes for this variability. 

4.6.1 On the possibility of a dynamically perturbed disc 

Considering that dynamically perturbed discs are found in other 
close binary systems, we have analyzed the possibility that the long 
cycle is produced by a tilted circumprimary disc whose normal axis 
precesses around the rotational axis of the gainer. In this case we 
should observe deeper CA when the disc edge hides the gainer, and 
shallower CA when a larger projected disc area is observed. How- 
ever, a tilted disc hypothesis presents severe problems: (i) it should 
produce large orbital changes in the equivalent width of the emis- 
sion lines proportional to the projected area and also large changes 
in the peak separation at high state, which are not observed, (ii) it 
is difficult to imagine how the optically thick disc might emit at 



Balmer lines and (iii) the hypothesis is not consistent with the sta- 
bility of the orbital LC solution (M12). In addition, the constancy 
of the V I R light curve at high and low state strongly argues against 
large scale dynamical perturbations including disc precession or 
disc wrapping. Beside the previous evidence, the conclusion that 
the light source of the long-cycle is not eclipsed (Ml 2) is not com- 
patible with the perturbed disc hypothesis. 

4.6.2 The cyclic wind and the long cycle 

We have found evidence that the long cycle in V393 Scorpii is 
produced by a wind injecting larger amounts of mass into the cir- 
cumbinary medium at high state. This conclusion is in agreement 
with the suggestion given by Ml 2 that circumbinary material is re- 
sponsible for the extra light during the high state. These authors 
give 9 x 10~ 9 Mq /yr for the mass transfer rate which means that 
6.2 x 10~ 9 Mq are injected into the interstellar medium every 
long cycle in steady state regime. This extra mass produces addi- 
tional line emission in low projected velocities and also additional 
continuum light by free-free emission. In this view it is the chang- 
ing wind that produces the long photometric cycle. The simultane- 
ous increase in Ha line emissivity near the donor and the gainer 
suggests a connection between the wind and the chromospherically 
active donor star. The wind could be caused by stream deflection at 
the disc hotspot. The possibility of magnetic fields connecting the 
donor and the disc and eventually powering the wind deserves to 
be studied, especially knowing that the binary is a transient radio 
source (Stewart et al. 1989). 

It is possible that the final cause for the long cycle is self- 
modulation of the mass transfer rate. As the gainer is rotating at 
critical rotation it cannot accumulate the extra mass received from 
the donor and this mass is finally expelled by the wind. Spec- 
troscopy shows that the high and low states in V 393 Sco corre- 
spond to changes in line and continuum emissivity, and probably 
density, of the wind, which could reflect changes in the mass trans- 
fer rate. Let's assume that the mass transfer rate is modulated by 
the photon flux arriving from the gainer. This flux is attenuated 
by intervening material when the wind is stronger, producing less 
irradiation flux and lower M^. This produces a decrease in wind 
strength and subsequently larger irradiation flux, bigger M2 and so 
on. This valve mechanisms could explain the long cycle, but theo- 
retical work is needed to demonstrate that it is a valid mechanism 
and stable on the long term, beside of taking into account the obser- 
vational evidence indicating a chromospherically active rather than 
an irradiated donor. 

The wind discovered in V 393 Sco might be driven by the 
same (still unknown) physical mechanism that the wind found in 
/3 Lyrae (Harmanec et al. 1996), but there are differences between 
f3 Lyrae and V 393 Sco that have been pointed out by Ml 2. These 
authors propose that f3 Lyrae is in an earlier evolution stage than 
V 393 Sco, and characterized by a larger mass transfer rate. 

4. 6. 3 On equatorial and polar outflows 

Our earlier claim for equatorial mass flows in V393Sco(M10) 
can now be tested with our new optical spectroscopic dataset. 
Our data suggest that line asymmetries, especially those seen 
around <E> = 0.5 and previously interpreted as equatorial mass 
flows through the L3 point, are part of a more complex orbital pat- 
tern. Although in this paper we have studied the asymmetries of 
the OI 7773 and He 1 5875 lines and M10 those of the He 1 1083 
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nm line, it is reasonable to assume that all these lines behave simi- 
larly during the orbital cycle. The new datasets show that the blue 
depressed wing remains during large part of the 2nd half cycle and 
cannot be interpreted in terms of an L3 outflow only (Fig. 22). 
The new pattern is best explained in terms of an anisotropic wind, 
emerging near the system barycenter in direction of the binary mo- 
tion. On the other hand, the interpretation of the red depressed wing 
near $ = 0.9 and 0.1 in terms of gas stream infall remains solid in 
this study, but the new observation of blue depressed wings around 
$ =0.15 is puzzling. A possible explanation for this feature is 
detection of material bouncing back from the stream/disc impact 
region. 

M10 favored equatorial mass flows as the cause for the long 
cycle in part due to the absence of long cycle variability in UV su- 
perionized lines. The resonance lines of Si IV, C IV and N V do not 
seem to be associated to the long-cycle phenomenon. It is then pos- 
sible that these lines are formed in the stellar wind emerging from 
polar regions of the B-type gainer (as revealed by line asymme- 
tries discussed by M10), but probably this polar-wind is not related 
to the disc-wind discussed in this paper, that is clearly involved in 
the long-cycle variability. The IUE observations are limited to only 
a few spectra, this unfortunate constraint impedes to have a clear 
picture of superheated regions during the long-cycle. 



5 CONCLUSIONS 

In this paper we have analyzed high resolution optical spec- 
troscopy of V 393 Sco with time resolution suitable for studying 
independently the orbital and long cycle variability. In general, 
the spectroscopic variability can be explained in terms of several 
circumstellar zones including a circumprimary disc, a gas stream, 
a hot spot and a modulated high latitude wind. We have shown that 
the long cycle is the manifestation of this variable wind, probably 
produced by the interaction of the gas stream and the optically 
thick disc surrounding the critically rotating B-type primary. Our 
main conclusions are: 

• We detected the donor in our spectra; its effective temperature 
was determined from a comparison with theoretical synthetic spec- 
tra, obtaining T2 = 7900 K. The donor RV was measured and the 
orbital elements determined; they are given in Table 5. We find a 
significant small orbital eccentricity for the system, but its meaning 
is not clear in an interacting binary with mass flows. 

• From the strength of diffuse interstellar bands we estimate 
E(B - V) = 0.15 ± 0.05, consistent with the M10 result. 

• We developed a method of spectral disentangling subtracting 
from every spectrum at every observing epoch the donor synthetic 
spectrum weighted by its fractional contribution given by our light 
curve model. This method turned to be adequate to remove the main 
spectral features of the donor and for studying the remaining com- 
ponents of the system. 

• At low stage the He 4471 and Mg4481 equivalent widths are 
consistent with line formation in an optically thick region with tem- 
perature 16.000 K. We argue that at low state these lines are less 
contaminated by emission and probe the pseudo-photosphere of the 
B-type gainer. 

• Residual spectra reveal several metallic emission lines from 
the donor, like Mg II 4481 and CI 6588. Doppler maps of these 
lines reveal a homogeneously-distributed emissivity around the 
secondary, not compatible with irradiation-induced emission but 
with a chromospherically active star. The highly variable nature 



of this emission along with the reported radio-transient nature of 
the binary rise the question on the probable magnetic nature of the 
donor. 

• Novel spectral features were detected: discrete absorption 
components (DACs), especially visible at blue-depressed 1 7773 
absorption wings during the second half-cycle, Balmer double 
emission with V/R curves showing "z-type" and "s-type" excur- 
sions around secondary and main eclipse, respectively, and H/3 
emission wings extending up to ± 2000 km s _1 . 

• DACs are a new and still enigmatic phenomenon; deeper stud- 
ies are needed to enlightening the nature of these features. We sug- 
gest that DACs could be Zeeman separated absorption lines formed 
in magnetized gas streams. This conjecture can be tested with spec- 
tropolarimetry. 

• Line asymmetries can be interpreted in terms of photon ab- 
sorption in the canonical gas stream, a large-scale wind and mate- 
rial bouncing back from the stream/disc impact region. 

• We notice that the relatively large widths of the Balmer and 
helium absorption cores are consistent with rotational broadening 
in absorbing material around a critically rotating gainer. 

• Balmer and He I emission are not eclipsed and in average the 
Ha violet peak is larger than the red one. At phases of maximum 
disc visibility and on high state the emission profiles are single- 
peaked. These findings suggest that Balmer and He I emission are 
produced in a high-latitude bipolar wind. 

• The extended Balmer emission wings can be the result of 
Thompson scattering in dense regions of this wind. These regions 
are in the binary hemisphere containing the stream-disc interaction 
region and could be formed by stream deflection at the hot spot. 

• As cause for the long cycle we exclude a dynamical instabil- 
ity like disc precession or disc wrapping. This conclusion is based 
on the stability of the V/R curve, the stability of the photometric 
orbital light curve, the variability of spectral features like peak sep- 
aration and most importantly the fact that the long cycle is produced 
by a non-eclipsed source. 

• We interpret the long-cycle in terms of a bipolar wind whose 
emissivity is modulated with the long period. This wind produces 
rather stationary Balmer, OI, Sill and He I broad emissions. On 
high state the wind strength is larger than on low state. We argue 
that an increase in wind strength is the cause for the larger line 
emission and brighter and redder continuum on the high state. 

• It is very likely that the wind is the result of the interaction 
between the gas stream and the optically thick and massive circum- 
primary disc. A donor magnetic field could play an important role 
in the energetic of the phenomenon and acting as a driving mecha- 
nism for the wind. 
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